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. !" l ™ duct j° n - Analytic solutions for the response of planetaiy lithospheres to volcanic loads 
dSSh^Hnn n^fn° i!? 0 * 1 ® fau J ng f™* . in ( cr cIastic P late thicknesses [e.g., 1-4]. Predictions of the 

Sni VOl ffl C ! ? ads ‘ based on the a PP»cation of Anderson's criteria for 

faulting [5] to the results of the models, do not agree well with observations [6,7]. Such models 

d ° p e ^ css . state in l° ad its clf, but only suggest a state of horizontal compressive 

• . ■ . u er jj bese tpcdpls have considered only the effect of an instantaneously emplaced 

d ° 4 °i addre f s ^ Gme evolution of stresses, nor do they consider the effect of a load 
^£2?H 0 i VS 'i ^ fini ? element approach allows us to assign elements to the load itself, and thus 
permits calculation of the stress state and stress history within the edifice. The effects of episodic 
load growth can also be treated. When these effects are included, models give much better 
agreement with observations. 

Method. We use the finite element code TECTON [8,9] to construct axisymmetric models of 
volcanoes resting on an elastic lithospheric plate overlying a viscoelastic asthenosphere. We have 
implemented time-dependent material properties in order to simulate incremental volcano growth. 
Theviscoelastic layer was taken to extend to a sufficient depth so that a rigid lower boundary has 
«n/i S1 ,k ni ^ ant in fiuence on the results. The code first calculates elastic deformations and stresses 
and then determines the time-dependent viscous deformations and stresses. Time in the model 
scales as the Maxwell time x m in the asthenosphere. 

lithncnh at 1 An a volcano 25 km in height a " d 200 km in radius on an elastic 

lithosphere 40 km thick (parameters approximately appropriate to Ascraeus Mons). The volcano 

Ioad aments applied at intervals of 1000 t*. Contours of maximum deviatoric 

“ Af fu Uy-grown edifice at the conclusion of flexure (r = 3000 r m ) are shown in Figure 1. 
Note that the maximum stress occurs in the lower mid-section of the edifice. We adopt the 
convention that tension is positive. H 

Discussion. Simple models of plate flexure due to an instantaneous load predict a state of 
horizontal compressive stress in the plate beneath the load [1,4] with an increasing magnitude of 
TOmpressmn from the middle of the plate upward, and a maximum at the surface. Models with a 
succession of loads emplaced from the bottom up display a different distribution of stress The 
maximum horizontal stress occurs in the lowest (first) load increment and decreases upward with 
each later load increment. Thus, thrust faults (predicted by simpler models) associated with such 
stresses would be expected to occur only early in the evolution (or not at all) and would be covered 
by later units which remain unfaulted. Such a distribution of stresses could also affect the 
ocations and dimensions of magma chambers. A zone of horizontal compression within the 
edifice would inhibit ascending magma from reaching the summit and thus might cause magma to 
pool beneath it. Since magma will propagate perpendicular to the direction of least compfessive 
stress, radial propagation of magma in sheet dikes or sills might also occur. The effects of adding 
fault slip alon g a bas^ detachment between edifice and lithosphere (as proposed in models of the 
flank tectonics of Hawaii [10]) was shown to achieve an even greater reduction in horizontal 
compressive stress m the edifice [1 1]. Such a structure may control the growth and evolution of 
tectonic features surrounding volcanoes on both Mars and Earth [11] 

Stresses calculated from planetary loading models have been used, in conjunction with 
Anderson s criteria for faulting [5], to predict types of faulting that should be observed [1 2 4] 
Typically, a zone is found surrounding the load with both the least compressive stress ct, and the 
most compressive stress a 3 horizontal (assuming that near the surface, a„ is small and 
compressive). Application of Anderson's criteria yields a prediction of strike-slip faulting Such a 
region surrounding a large load has not heretofore been observed [6.7]. Given the complete stress 
field m an element in this region from the TECTON solutions, we can plot Mohr diagrams as in 
Figure 2, and compare them with a Mohr-Coulomb failure envelope. The Mohr circles grow in 


MODELLING TECTONICS OF LARGE VOLCANOES: McGovern PJ. and Solomon S.C. 


size as flexure proceeds, until the failure envelope is exceeded at about 70 Xm after first loading. 
The largest circle shows stresses at the conclusion of flexure. This stress state corresponds to that 
given by analytic plate models (in which flexure occurs instantaneously). Simply interpreting this 
final stress state is seen to be misleading for two reasons. First, Oi (= On) is much larger than the 
tensional strength of any rock. Use of a shear failure criterion (such as Mohr-Coulomb) for such a 
stress state is improper [7]. Use of a tensile failure criterion will predict the formation of tension 
joints which may be precursors of graben [7]. Second, consideration of the stress history indicates 
that at the time of first failure, the shear failure criterion is satisfied, with principal stress directions 
that predict circumferential normal faulting. Formation of these faults would relieve stresses and 
provide planes of low resistance along which further faulting with a similar orientation would be 
expected to develop. Thus, the "prediction" of a zone of strike-slip faulting surrounding a 
lithospheric load, based on simple instantaneously-loaded plate flexure models, is m error. 

Conclusions. Simple interpretations of the results from instantaneous load models lead to 
predictions of lithospheric stress fields that do not match observations. Consideration of the 
history of stresses in the lithosphere and within the load itself helps to resolve these discrepancies. 
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Figure 1. Contours of maximum deviatoric stress (MPa) in the volcanic edifice described above. 



Figure 2. Mohr diagram for element at the top of the lithosphere (stresses calculated at element 
center: r = 265 km, z = -2.5 km). Diagonal line is Mohr-Coulomb failure envelope for basalt [12] 
with angle of internal friction 4 > = 49° and cohesion c = 38 MPa. In order of increasing radius, the 
circles represent stress state at t = 30, 70, 1000, 2000, and 3000 x m . When failure envelope is 
first exceeded (r = 70 x m ), Oj = and C 3 =Gzz- At final time shown, cj; - G„ and G 3 = Oqq. 


